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the occupation probability P{k) with which an urn has k balls by using the concept 
of statistical physics of disordered systems. In the disordered urn model belonging 
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We also discuss an interpretation of our results for explaining of macro-economy, in 
particular, emergence of wealth differentials. 
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1. Introduction 

A lot of techniques and concepts of statistical mechanics of disordered spin systems, in 
particular, the replica method originally used to analyze the thermodynamics of spin 
glass model by Sherrington and Kirkpatrick [1], have been applied to various research 
fields beyond conventional physics, i.e. information processing [2], game theory [3] and 
so on. The exactly solvable mathematical model, which describes these problems, is 
categorized as mean-field class [1]. 

On the other hand, as another exactly tractable model, in 1907, Paul and Tatiana 
Ehrenfest published a paper corroborating Boltzmann's view of thermodynamics [5]. 
Their urn model has been defined by Kac ^ as an exactly solvable example in statistical 
physics. While it has also been criticized as a marvelous exercise too far removed from 
reality, their urn model has been applied to modern problems such as complex networks 
O [8] or econophysics P |10], etc. For instance, based on extensive simulations of 
the Lennard- Jones fluid requiring in part a parallel computer in Juelich, an Italian- 
German team has shown that the prediction of the Ehrenfest urn effectively describes 
the behavior of the gas phase [H]. Moreover, it has been revealed that the mathematical 
structure of equilibrium state of the urn model [12] is similar to the zero-range process, 
which has been widely investigated in research fields of non-equilibrium statistical 
physics [13]. 

Recently, in the research field of complex networks ^ ^Sj, Ohkubo et al. [H] 
proposed a network model based on 'Ehrenfest class urn model' to explain how the 
complex network gets scale-free-like properties, where 'Ehrenfest class' means that each 
urn has distinguishable balls. In the model, each urn corresponds to a node in graph 
(network) and the number of distinguishable balls, k, in each urn is regarded as degree of 
nodes. For this model system, they succeeded in deriving the scale-free-like properties 
~ fc~^(log/c)~^ in the probability of the degree of nodes by the usage of the replica 
symmetric theory [15]. In addition, the similarity between the disordered urn model and 
the random field Ising model [16], and the condensation phenomena in the disordered 
urn model have been investigated |17j . 

We here note that there are a lot of works in which the power-law behavior and the 
condensation phenomena in urn models have been studied [121 [18]. For example, in the 
zeta urn model [12] , the power-law behavior in the probability of the number of balls, 
i.e., the occupation distribution, stems from a power-law form of the Boltzmann weight. 
However, when we attempt to describe various problems in the real world, we should 
take into account the disorder and treat the urns as a heterogeneous system. Mainly, 
the previous models which cause power-law behavior in the occupation distribution do 
not contain any disorder, and hence it would be important to investigate 'disordered' 
urn models which cause the power-law distribution. 

In this paper, we propose two disordered urn models in which quenched randomness 
is important for generating the power-law behavior. One of them belongs to the 
Ehrenfest class, and the other corresponds to the Monkey class, in which each urn has 
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indistinguishable balls. In particular, for the Monkey class urn model, we investigate 
a real-space condensation phenomenon in which macroscopic number of balls are 
condensed into only one urn. The occupation probability P{k) with which an urn 
has k balls is calculated analytically, and furthermore, the critical density pc for a given 
temperature is evaluated. As the result, it is shown that the occupation distribution 
function P{k) changes its scaling behavior from the exponential e~'^-law to the 

power-law in large k regime. 
This paper is organized as follows. In the next section [2], we introduce the general 
formalism for the urn model with an arbitrary energy function. Although there are 
several analytical treatments for disordered urn models fi5[ IT9]. we give the formalism 
for the disordered urn model with an arbitrary energy function with a different point of 
view, and additionally, the analytical treatment makes this paper self-contained. This 
analytical treatment contains both Ehrenfest and Monkey classes as its special cases. 
We explain the relation between the saddle point, that determines the thermodynamic 
properties of the system, and the chemical potential. With the assistance of this general 
formalism, we provide an analysis for a special choice of the energy function, which is 
an example of the Ehrenfest class in section [3l We discuss the condition on which the 
power-law appears in the tail of the occupation probability for the model. In section 
m we show that the condensation occurs for the special case of the Monkey class with 
disorder, and a heavy tailed power-law emerges in the occupation probability. In section 
m we provide a possible link between our results and macro economy, in particular, 
wealth differentials. Last section is a summary. 

2. General formalism for urn models with disorder 

Let us prepare urns and M balls (M = pN) and consider the situation in which the 
A^ urns share the M balls. Then, we start our argument from the Ehrenfest class urn 
model [12] in which each ball in urns is distinguishable. For the mathematical model 
categorized in the Ehrenfest class, the Boltzmann weight p{ni) that i-th urn possesses 
rii balls is given by 



where E{ei, ni) is an energy function, a disorder parameter for urn i, and (3 the inverse 
temperature of the system. The factorial n\ stems from the property of the Ehrenfest 
class [12] . The only point of our analysis which is different from [12] is quenched disorder 
ei appearing in the energy function. For the Ehrenfest class, the probability that an urn 
specified by the disorder parameter ei possesses k balls is given by 



pijii) 



exp [-l3E{ei,ni)\ 



(1) 




Hi — M — l \ log z 
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and used the Fourier transform of the Kronecker-delta: 
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(3) 
(4) 



to introduce the conservation of the total balls : ni + - ■ ■+nAr = M = pA^ into the system. 
In order to calculate -P(A;), we take the configuration average of /fc''^^({ei/i})- While one 
can calculate the configuration average by means of the replica method [15], it has been 
revealed that the mathematical structure of the disordered urn model is related to that 
of a random field Ising model [16j. Hence, we here use the law of large numbers and 
simplify the calculation. To calculate the average of the quantity exp[log /^^^^({cj/i})] 
over the configuration {ei, ■ ■ ■ , eTv}, we consider the Tayler-expansion: 

(exp[log/i^^)({e,/i})])|i,2,...,iV} = 1 + (log/fH{e,/i})){i,2,.,,v} 



+ 2((l°g/^^({^^/l}))')av■,^} + ■■ 



(5) 



where (■ ■ ■){i,2, --,Af} means the configuration average over {ei, ■ ■ ■ , e^^. Here we assume 
that the observable log/^^^^({ei/i}) for a given realization of configuration {ei, ■ ■ ■ , eAr} is 
almost identical to the average (log/^'^^''({ej/i})){i^2,- -,Af}- In other words, the deviation 



is vanishingly small as ((log /i''^({ei/i}))2)|i 2,...,7V} - (log /fc''^({ei/i}))|i 2,-,iv} - in the 
thermodynamics limit. By using the assumption, we rewrite (j5]) as follows. 

(exp[log/i-)({e./i})]){i,2,.,,v} ^ 1 + (log/fH{e,,/i})){i,2....,iv} + ^(log/^^^({e./i}))?i,2,...,iv} 

+ ■■■ = exp[(log/^^({6,/i}))|i,2,...,TV}] (6) 
This replacement of the configuration average reads 

(/i"^({e./i})){i,2,...,TV} = (exp (log/i"H{evi})))^, , ...^^j ^ exp ((log /(-^({e,/ J))|,,2,...,tv}) 

= exp (^(log[e-^^(^i''=V^!] + (log Z2){2,3,.,7V} - (log Zi){2,3,...,7V})^^ J . 

In the thermodynamic limit 00, (log Z2){2.3,...,Ar} is evaluated as 



(7) 
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dr ^ 
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(A;_piV-l)log^ + ^ log^: 
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ni=0 
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where we used the law of large numbers, and in the final line the saddle point method 
was used. (■ ■ ■) means the average over only e, namely, (■■•)= /(■ ■ ■)D{e)de. In the next 
two sections, we consider specific choices of D{e) to evaluate the occupation probability 
distribution concretely. Using the same way as the Z2, Zi, which is obtained by the 
normalization condition of fk^\ namely, J2kLofk^^ — ^ with equation ([2]), is rewritten 
as 

(log ^l){2,3, -,7V} 

= log E 4 - (PA^ + 1) log ^.2 + (iV - 1) ( log E )• (9) 

fc=0 ^ n=0 ^' / 

We easily find z^g = Z2s because the first terms for each saddle point equation ([8]) or ([9]) 
are vanishingly smaller than the other two terms in the limit of N —>■ 00. 

Thus, we obtain the saddle point equation with respect to Zs = Zsi = Zs2 and the 
occupation probability P{k) = (f^^) that an arbitrary urn with the energy function E 
at inverse temperature jS has k balls are given by 

-1 /^n=(}"' (n-iy. \ , , 



n 



=0 



and 



e-'3-E(.,fc)_^fc 



Z^n=0 n! 

respectively. It should be noticed that the above saddle point equation for the Ehrenfest 
class urn model (llUp is now rewritten in terms of chemical potential 

fx = p-HogZs (12) 

as 

g — /3[E(e,n)— n^i] 



Then, we define the probability pn that an arbitrary Ehrenfest class urn with energy E 
has n balls by 

0w(e,n) e-/3[i;(..n)-nM] 

From this formula of the probability Pn with the effective Boltzmann factor 4'E,fi,i3{^, n), 
the equation (IT^ means that the ratio M/N corresponds to the average number of 
balls put in an arbitrary urn: p = J2'i^=o^4>E,iJ,,i3{^,n)/ J2^=o(pE,iJ.,/3{^,n), and its value is 
controlled by the chemical potential fi through the equation f[T51) . Then, the chemical 
potential /x and the saddle point Zg are related through the equation (fT2|) . Therefore, 
when we construct the system so as to have a density p, the corresponding saddle point 
Zs is given by f|T3l) . As the result, the chemical potential p that gives p is determined 
by the relation (|T2|) . 
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Thus, our problem is now to solve the saddle point equation 

/ E^=o^0w(e,^) \ ^^^^ 



and to calculate the following averaged occupation probability for the solution Zg = e^^ 
of the equation f|T5|) : 



Now it is time for us to stress that the Ehrenfest or Monkey class is recovered if we 
choose the effective Boltzmann factor 0£,/^,/3(e, n) as follows [T2] . 

, , J exp [— /5(i?(e, ra) — n/i)] (Ehrenfest class), , , 

0E,M,/3(e,ri) = j [-/5(E(e,n) - n/i)] (Monkey class). ^ 

It should be noted that in our formalism, the distinction between two models only comes 
from the difference of the effective Boltzmann factor (ITTl) . We here also comment on 
the effect of the disorder to the dynamics of urn models. In the uniform case without 
disorder, all urns are equivalent, and we do not distinguish each urn in the Monkey 
class. On the other hand, in the disordered urn model, each urn i has own label ej, 
and hence they are distinguishable at least in principle. However, we assume that this 
heterogeneity of urns does not change the "Monkey" nature in the Monkey class; we 
cannot see the disorder parameter assigned to each urn from outside, and the dynamics 
(so called 'box-to-box choice') is not changed. 



3. Ehrenfest class urn model with disorders 



As an demonstration of the Ehrenfest urn model whose thermodynamic properties are 
specified by equations (ITOl) and (fTTi) . we introduce a kind of disordered Ehrenfest urn 
models and consider the condition on which the power-law appears. To this end, we 
choose the energy function E{e,n) as 

E{e,n) = - en, (18) 

where e means an urn-dependent disorder of the system, and we here assume that e 
takes a value in the range [0, 1] randomly, that is, D{e) = 6(e) — 9(e — 1) with the 
step function 0(- ■ ■). Of course, we might choose the other distribution -D(e), however, 
the main issue of this section is to examine whether the power-law appears in the 
occupation probability distribution when we introduce the disorder, and for this purpose, 
the simplest choice is enough. The tendency of this energy function to force each urn of 
the system to gather balls as much as possible results in the fact that the rich get richer 
as its collective behavior. We should mention that in |19], the co-called Backgammon 
model [20] described by the cost function E{ei,ni) = —etSmfi for each urn was studied. 
In the model, the cost decreases if and only if each urn is empty. In this sense, the 
model we deal with in this section is regarded as an opposite situation of our model. 
Therefore, it is interesting to investigate whether there exists any significant difference 
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or the similarity between two models from the view point of the occupation probability 
distribution. The extensive studies concerning this issue will be our future studies. 

For this choice of the energy function ( fT8|) . the saddle point equation ( ITOl) leads to 



' e/3 - 1 

From equation fill I) , the occupation probability for the choice (1181) . P{k), is given by 

zk ,1 



(19) 



Pik) 



k\ Jo 



deexp [Pek — Zge 



(20) 



In following, we evaluate the above occupation probability. We first show the phase 
diagram that indicates the area of the power- law behavior. In Figured] (left), we show 




= 2 ■■ *u 

= 5 ; *:„ 

= 2 ■■ *d 
--3 :k„ 
= 5 ; 




Figure 1. The p-k phase diagrams concerning the heavy tailedness of the occupation 
probabihty. In the phase diagram for (3 — 3 (left panel), the occupation probability 
follows non-power, exponential-law. The power-law regime exists for p > ph = 
2.697423. In the right panel, we plot the phase boundaries kd,ku as a function of 
density p for several values of inverse temperature, namely, /3 = 2, 3 and 5. 



the phase diagram for the case of /? = 3. The shaded area in this figure means non- 
power, exponential-law region. From this figure, we find that while the upper bound 
(the cut-off) fcu and the lower bound k^ increase as the density p increases, the cut-off 
fcu increases much more quickly than the lower bound k^- As the result, the heavy tailed 
power-law fc~^-region is broadened by increase of the density p. In the right panel of 
Figure [1], we display the inverse temperature dependence of the area of the power-law. 
As temperature increases, the area of the power-law shrinks to zero. Then, we have 
a Poisson-law in whole region of the phase diagram in the high-temperature limit of 
(3 = 0. 

We explain the detail of the evaluations as follows. We first consider the high- 
temperature limit P = 0. For this case, the saddle point ( |T9|) leads to Zg = p and we 
obtain 

(21) 



Pik) = -j^ exp (-p) , 
which is nothing but a Poisson distribution. 
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For finite temperature f3 > 0, the occupation probability P{k) is rewritten as 



P{k) 



1 



e/^-l 



oo 



, oo 



(22) 



where we defined the incomplete Gamma function by Tj^{a,b) = jj^t^ ^e ^dt. In Figure 
121 we plot the occupation probability P{k) for several values of p for /5 = 3, Zs = 0.5, 1, 2 
and 3, namely, p = 3.180923,6.361846,12.723691 and 19.085537. From this figure, we 





Figure 2. The Log-Log plot of occupation probability P{k) as a function of 
k for several values of z^, namely, Zg — 0.5,1,2 and Zg ~ 3, namely, p = 
3.180923,6.361846, 12.723691 and 19.085537. We set the inverse temperature /3 = 3.0. 
The right panel shows the behavior of <i>(/c). 



find that P{k) obeys a power-law ~ k~^ in the intermediate regime of k and there exists 
a cut-off value from which the distribution decays exponentially. To see the existence 
of the cut-off value explicitly, we rewrite the distribution P{k) as 



0, 



(3p 



f3peP 



oo 



(23) 



e/3p - 1^ 

In Figure [2] (right), we plot the function $(A;) for several values of Zs a.t (3 = 3. 
We easily find that in the range of [k^, ^u] defined in terms of the function $(A;) as 
kd = min{/c|$(A;) = 0}, fcu = max{A;|$(A;) = 0}, the occupation probability P{k) follows 
a power-law distribution ~ k~^. Thus, we specified the region (p, k) in which a power- 
law heavy tail appears in the occupation distribution as shown in Figure [H 

As shown in this section, for the Ehrenfest class disordered urn model, we could 
clarify the control parameters of the system for which the heavy tail power-law emerges 
in the occupation probability. In the next section, we consider a Monkey class urn model 
with disorders. 



4. Bose-Einstein condensation and emergence of the heavy tail 

In the previous section, we evaluated the asymptotic form of the occupation probability 
P{k) for the urn model of the Ehrenfest class with energy function E{e,n) = —en. 
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Obviously, from the view point of the energy cost, it is a suitable strategy for each 
urn to gather balls as much as possible. In that sense, this case should be referred by 
the concept the rich get richer in the context of social networks. However, by using 
the general definition of the problem, we freely choose the energy function for both the 
Ehrenfest and Monkey classes. 

In this section, for the Monkey class urn model whose thermodynamic properties 
are defined by equations f|T5|) . f|T6|) and f|T7j) . we evaluate P{k) for a specific choice of 
energy function E{e,n). We choose the energy E{e,n) as 

E{e,n) = en (e > 0). (24) 

We should notice that for this simple choice of the energy function, the urn labeled by 
e 7^ is hard to gather the balls. On the other hand, the urn with e = energy level 
easily gathers the balls. The urn model having this type of energy function does not 
agree with the concept the rich get richer. Nevertheless, we use the energy function ( 12^ 
because as we shall see below, a kind of condensation with respect to the urns occurs 
for this choice of energy function, and as the result, the power-law in the tail of the 
occupation probability emerges. 

For a given choice of -D(e) as the density of state, namely, degeneracy of the energy 
level of the urn, we rewrite the saddle point equation ( |T5l) as follows: 

To proceed to the next stage of the calculation, we choose the density of state D{e) 
explicitly as 

D{e) = eoV~e, (26) 

where is a constant. Although we chose the above form, more general setup of the 
argument by choosing -D(e) = eoe" is possible. We shall discuss the result for this kind 
of generalization later on. Then, the equation fl25|) is rewritten by 

p eo^/ede 

P= L .-ie/^--l +^-°' (27) 
where Pe=o means the density of balls in the urn labeled by the zero-energy level e = 0. 
We should notice that the second term appearing in the right hand side of equation 
fl27|) . namely, Pe=o vanishes in the thermodynamic limit N —>■ oo when a condensation 
does not occur. In other words, when the condensation arises, the second term, Pe=o, 
becomes from zero to a finite value; this means that an urn with e = becomes to have 
a macroscopic number of balls. 

In following, we show the system undergoes a condensation and investigate the 
behavior of the system when the density p increases beyond the critical point pc for a 
given finite inverse-temperature /3. 

• Before condensation: p < Pc 

By a simple transformation jSe = x, the equation fl271) is rewritten in terms of the 
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so-called Appeli function (see e.g. [H]) bn{zs) as follows: 



P = ^r^/^&3/2(..), (28) 

where the Appeli function (see e.g. [21]) bn{zs) is defined by means of the Gamma 
function Tin) as 



b (z)-^ r^^^ (29) 
~ T{n) Jo ^ 

We should keep in mind that &3/2(-2s) < ^3/2(1) = C(3/2) = 2.6... is satisfied 
(6„(1) = Y^kLi = ({n)). The solution of the saddle point equation fl29l) possesses 
a solution Zs < I. 
At the critical point: p = Pc 

The critical point at which the condensation occurs is determined by the radius of 
convergence for the following partition function: 

00 00 

Z = J2 -^re"^'" = ^(e-^^+i°s^=)", (30) 

n=0 n=0 

namely, Zg = 1 for e = gives the critical point. Therefore, substituting Zg = I for 
a given density level p, the critical point pc, above which the condensation occurs, 
is obtained by 



Pc = ^/?-/^ 63/2(1). (31) 

• After condensation: p > Pc 

For p > Pc, the saddle equation of (129|) no longer has any solution as 2;^ < 1. 
Obviously, for the solution > 1, the partition function diverges. Then, we should 
bear in mind that the term p^=o in fl27|) . which was omitted before the condensation, 
becomes C(l) object and the saddle point equation we should deal with is not ( l29ll 
but fl27j) . As the result, the equation (1271) has a solution Zg = I even for p > Pc and 
the number of balls fc* in the condensation state increases linearly in p as 

K = N{p-pc), (32) 

whereas the number of balls in excited states reaches k = Npc- 

Thus, we obtained the saddle point Zg for a given density and inverse-temperature. We 
found that the condensation is specified by the solution Zg = 1. 

We next investigate the density dependence of the occupation probability through 
the saddle point. For the solution of the saddle point equation Zg, the occupation 
probability at inverse temperature (3 is evaluated as follows. 

-<^-' ^ f (eS^) ^ T^-^'^ - - 

The above occupation probability is valid for an arbitrary integer value of k for k > 1. 
We plot the behavior of the occupation probability P{k) in finite k regime in Figure El 
In this plot, we set Eq = 1 and Zg as Zg = 0.1,0.8,1.0, and the inverse temperate 



Power-law behavior and condensation phenomena in disordered urn models 



11 




Figure 3. The behavior of the occupation probabihty ((33|) in non-asymptotic regime. 
We set eo = 1 and Zs as Zg = 0.1,0.8,1.0, and the inverse temperate is /3 = 1. The 
inset of the figure shows the asymptotic behavior of the occupation probability 
as Log-Log plots for the cases of — 0.1, 0.5, 0.7, 0.99 and 1. 



is /? = 1. In the inset of the same figure, we also show the same data in Log- 
Log scale for the asymptotic behavior of the probability P(/c) for several values of 
Zs, namely, = 0.1,0.5,0.7,0.99 and 1. From this figure, we find that the power- 
law emerges when the condensation is taken place for p > p^. The numerical 
analysis of the occupation probability (l33l) in the limit of A; — oo is easily confirmed by 
asymptotic analysis of equation (133|) . We easily find that the asymptotic form of the 
wealth distribution P{k) behaves as 

P{k) = r^/^z'^^il - Zs)eoT{3/2) k'^^^ + ^(3-^/h''^+^eoT{3/2) k-'''^ + 0{k-'"^). (34) 

We also should notice that a macroscopic number of balls k^, is gathered to a specific 
urn with energy level e = when the condensation occurs. As the result, the term such 
as ~ {1/N)5{k — k^) should be added to the occupation probability. 
Let us summarize the results as follows: 

p,,, _ / 2^A:-3/2exp[-Hog(l/^.)] {p < p. : z, < 1), 

- I 3^ ^_5/. + ^^(^ _ (p > : = 1). ^^^^ 

The scenario of the condensation is as follows. For a given p < pc, one can find 
a solution Zg < 1 for the saddle point equation (1271) . and hence the second term 
of equation (!27j) is zero in the thermodynamic limit. Then, for non-condensed Np 
balls, the occupation probability follows ~ k~^^'^ e~'^-law. Namely, urns possessing a 
large number of balls do not appear due to the repulsive force as E = en. When 
p > Pc, the saddle point Zg is fixed as = 1; if > 1; the first term of the saddle 
point equation flTTI) has a singularity. Therefore, in order to avoid the singularity, the 
second term of the saddle point equation (127|) becomes from zero to a finite value. As 
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the result, the occupation probability is described by the /c~^/^-law with a delta peak 
which corresponds to an urn of e = gathering the condensed N{p — p^) balls. This 
corresponds to the condensation phenomena in the disordered urn model. In particular, 
the occurrence of the condensation in the disordered urn model treated in the present 
paper is characterized by the transition from the exponential-law to the heavy tailed 
power-law. We also mention the effect of disorder on the power-law behavior of the 
occupation probability. We easily find that the power-law behavior disappears when 
one cancels the disorder of the system by choosing the density of the energy such as 
D{e) = 6{e — e) (e is a constant). This fact means that the disorder appealing in the 
system possesses a central role to make the occupation probability to have a power-law 
behavior. 

We should notice that in the above argument, the solution Zg = I that indicates 
the condensation does not change even if we choose the density as D{e) = eoe", (a > 0). 
For this choice, the pc is given by 

Pc = eof3-'-- (36) 
JO e^ — 1 

Then, one obtains the following occupation probability 

P(k) = '■"'It'' - ''"'f i? ^ " + 1)--'- (37) 



At the end of this section, we should mention the result for the uniform distribution of 
e, that is, the case of a = leading to D{e) = 0(e) — 0(e — 1). For this choice, we have 
the following occupation probability 

k fc+1 

P(,).i|(l-e-»)-^(l-e-«-"). (38) 

Then, beyond the critical density Pc = Jode/{e'^' - 1) = /3~^En=ii^ - e"^")/n, the 
occupation probability P{k) behaves as 

Therefore, after the condensation, the crossover from the fc~^e~^'^-law to the /c~^-law is 
observed around k ~ l3~^ and as the result, the power-law heavy tail appears. 



5. Interpretation from a view point of macro economics 

In this section, we reconsider the results obtained in the previous sections from a view 
point of macro economics. It is easy for us to regard the occupation probability as wealth 
distribution when we notice the relations: balls - money and urns - people in a society. 
In following, we attempt to find an interpretation of the condensation and the emergence 
of the Pareto-law [22] in terms of wealth differentials [23l[22[25l[26l[271[28l[29l[30l[31]. 

In section [3l we devoted our analysis to extremely large income regimes (the tail 
of the wealth distribution), however, it is quite important for us to consider the whole 
range of the wealth. As reported in |30], the wealth distribution for small income regime 
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follows the Gibbs/Log-normal law and a kind of transition to the Pareto-law phase is 
observed. For the whole range distribution of the wealth, the so-called Lorentz curve 
[321 [33l is obtained. The Lorentz curve is given in terms of the relation between 
the cumulative distribution of wealth and the fraction of the total wealth. Then, the 
so-called Gini index [32l [331 EH [35] , which is a traditional, popular and one of the most 
basic measures for wealth differentials, could be calculated. The index could be changed 
from (no differentials) to 1 (the largest differentials). For the energy function flMl) in 
the previous section, we derived the wealth distribution for the whole range of incomes 
k. In this section, we evaluate the Gini index analytically. 

As we mentioned above, the Lorentz curve is determined by the relation between 
the cumulative distribution of wealth X{t) = //^^.^^ P{k)dk and the fraction of the total 
wealth Y{t) = . kP{k)dk/ kP{k)dk for a given wealth distribution P{k). For 
instance, the Lorentz curve for the exponential distribution P{k) = 'ye~'^^ is given by 

F = X + (l-X)log(l-X). (40) 

We should notice that the Lorentz curve for the exponential distribution is independent 
of 7. 

For the power-law distribution P{k) = (7 — 1) A;"''' (7 > 1), we have 

Y = l-{1-X)^. (41) 

as the Lorentz curve. This curve depends on the exponent 7. In Figure [H we plot the 
Lorentz curve for the exponential distribution (HOl) and the power-law distribution (HTi) 
with several values of 7. 




0.2 0.4 Y 0.6 0.8 1 0.2 0.4 Y 0-6 0.8 1 



Figure 4. The left panel: The Gini index is obtained as an area between the perfect 
equality line Y = X and the Lorentz curve (e.g. Y = X + {1 — X)log(l — X) 
for exponential distributions). The right panel shows the Lorentz curves for the 
exponential distribution (|40|) and the power-law distributing ([4T|) with several values 
of 7 = 2.5 and 3. 



Then, as shown in the left panel of Figure [U the Gini index G is defined as an area 
between the perfect equality line Y = X and the Lorentz curve. This quantity explicitly 
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reads 

/■I roo fix 

G = 2 (X - Y)dX = 2 / (X(t) - Y{t)) ■ ^dt (42) 

and we have G = 1/2 |33l [34] for the exponential distribution and G = 1/(27 ~ 3) for 
the power-law distribution. As the occupation probability distribution fl33l) is defined 
for k > 1, one can evaluate the Gini index as a function of the saddle point Zg- In Figure 
El we plot the Lorentz curve (left) for several values of Zg. In the right panel, the Gini 
index G{zs) is shown. We find that the index approaches to 1/2 as Zs 1. 

From the argument in the previous section, we easily find that the occupation 
distribution for Np^ non-condensed balls beyond the critical point is modified such as 
~ k~^"~^'^^ by choosing the density of the energy D(e) = eoe". Namely, for the Pareto- 
law distribution ~ k~^°'~^'^\ the Gini index leads to G = l/(2a + 1). Therefore, the 
condensation is specified by the change of the Gini index from G = 1/2 to l/(2a -|- 1). 
However, we should keep in mind that the Gini index itself has less information about the 
differentials than the wealth distribution. For example, the Gini index for a = 1/2 of the 
Pareto-law P{k) ~ /c^^^'^^) gives the same Gini index as the exponential distribution. 
This fact stems from the definition of the Gini index G, that is, G is defined as an 
area between Y = X and the Lorentz curve. It could be possible to draw lots of the 
Lorentz curves that give the same area (the same Gini index). As we explained above, 
it should be noted that actually the Gini index is one of the measure for the earning 
differentials, however, the wealth distribution is much more informative than the Gini 
index. Although the Gini index is less informative than the distribution, however, for 




0.2 0.4 X "-^ "-^ ^ " 0-2 f-4 2, "-6 1-8 1 



Figure 5. The left panel: the Lorentz curve for ([55)1 . The right panel shows the Gini 
index for several values of Zs ■ 

a given real (empirical) data xi < X2 < ■ ■ ■ < xn {x denotes the amount of money 
for example. Such empirical data are not always massive enough for us to specify the 
distribution), one can evaluate it as a statistics by G = (l/X^yu) Z]ili(2^ — N — l)xi 
with /iiV = Y.^=i p5] . Therefore, it is helpful for us to use the Gini index to compare 
the earning differentials between different countries (the population of each country 
should be different and of course it is finite N < 00). In our analysis of this paper, the 
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distribution was analytically obtained in the tliermo-dynamic limit because we treated 
an ideal CclSG clS db society. Nevertheless, even if we encounter more realistic situation for 
which the analytical evaluation of wealth distribution is very tough, one can evaluate the 
earning differentials via the Gini index by computer simulations for finite population A^. 
Then, one can investigate the earning differentials by comparing the numerical results 
with the analytical expressions obtained in this paper. 

6. Summary 

In this paper, we investigated equilibrium properties of disordered urn model and discuss 
the condition on which the heavy tailed power-law appears in the occupation probability 
by using statistical physics of disordered spin systems. We applied our formalism to two 
urn models of both the Ehrenfest and Monkey classes. In particular, for the choice of 
the energy function as E{e,n) = en with density of state -D(e) = eoe"" for the Monkey 
class urn model, we found that above the critical density p > pc for a temperature, 
the condensation phenomenon is taken place, and most of the balls falls in an urn with 
the lowest energy level. As the result, the occupation probability changes its scaling 
behavior from the exponential k~^°'~^^^ e~'^-law to the k~^°'^'^^ power-law in large k regime. 
We also provided a possible link between our results and macro economy, in particular, 
wealth differentials. 

Of course, there might exist the other urn models showing the power-law behavior 
after the condensation. In fact, we find such a case in a recent study on the Ehrenfest 
urn model [I7j, in which the occupation probability follows a Poisson-law when the 
condensation occurs. Although we provided a piece of evidence to show that the 
power-law behavior in the occupation probability distribution takes place after the 
condensation for several restricted cases of the cost function, it is not yet clear whether 
the condensation always causes the power-law or not. The nature of the link between 
them will be a central problem to be clarified in future. Therefore, as one of our future 
studies, it might be important to investigate the universality class of urn models that 
shows the power-law behavior in the occupation probability beyond the critical point. 

We hope that various versions and extensions of the disordered urn model, including 
Backgammon model [201 IS] 5 could be applied to research area beyond conventional 
statistical physics. 
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